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Abstract 

Due to inter-subunit communication, multisubunit enzymes usually hydrolyze ATP in a 



On ' concerted fashion. However, so far the principle of this process remains poorly understood. In 

OO ' 

rn 

O ' this model, we assume that the binding of ATP will change the potential of the corresponding 



this study, from the viewpoint of statistical thermodynamics, a simple model is presented. In 



enzyme subunit, and the degree of this change depends on the state of its adjacent subunits. 
The probability of enzyme in a given state satisfies the Boltzmann's distribution. Although it 



j^ I looks much simple, this model can fit the recent experimental data of chaperonin TRiC/CCT 

well. From this model, the dominant state of TRiC/CCT can be obtained. This study provided 
a new way to understand biophysical processes by statistical thermodynamics analysis. 



Introduction 

In living cells, many enzymes consist of several subunits, each of which can bind and hydrolyze 
ATP by itself.-"^ Because of the positive or negative cooperativity among the enzyme subunits, 
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this ATP hydrolysis process is usually studied by Monod-Wyman-Changeux (MWC) model or 
Koshland-Nemethy-Filmer (KNF) model. ^"^ However, recent experimental measurements of the 
mammalian type II chaperonin TRiC/CCT, which consists of two ring-shaped cavities with lids 
composed of eight different subunits each, found that the bound ADP number presented on a 
single enzyme cannot be well described by the usual cooperative models, MWC model or KNF 
model.- 

In experiments,- the stoichiometry of hydrolyzed ATP (in the form of either bound ADP or 
the transition state mimic, ADP-AlFx) in a single TRiC enzyme is firstly measured, and then the 
hydrolyzed ATP number distribution is extracted. In this study, we assumed that all the bound ATP 
molecules will be hydrolyzed, and so the number distribution of hydrolyzed ATP is just the number 
distribution of ATP bound to one enzyme. As had been concluded,^ many standard cooperative 
models can be ruled out by the measured shape of ATP number distribution. In this study, we will 
not try to give a modified and so more complicated MWC or KNF model to fit the experimental 
data. On the contrary, we will present a simple thermodynamic model which is only based on 
Boltzmann's law. 

Results and discussion 

In our model, the enzyme TRiC/CCT is regarded as a double-ring with its two rings stack 
perfectly on one another (see [I]^~— ). Each ring consists of 8 CCT subunits which can bind and 
hydrolyze ATP by itself. The potential of each empty subunit (i.e. without ATP molecule) is 
denoted by e, and the potential of each free ATP molecule is denoted by ji- The (information) 
entropy of ^(> 1) free ATP molecules is denoted by 7oln^. For the ATP bound subunits, the po- 
tential of the ATP-subunit complex depends on the state of its adjacent subunits. For convenience, 
the two subunits which lie in different rings but be adjacent to each other, i.e. the two subunits 
with the same label number (see the schematic depiction in [T), are said to be mirror subunit to 
each other. In the model, we assume that the two rings of TRiC are perfectly stacked, i.e. each 



subunit has and only has one adjacent subunit which lies in the different ring. For one ATP bound 
subunit, if its mirror subunit is empty and there are total / ATP bound adjacent subunits, the poten- 
tial of this ATP-subunit complex is denoted by a,-, for < / < 2. On the other hand, if its mirror 
subunit is ATP bound, and there are altogether i ATP bound adjacent subunits, the potential of this 
ATP-subunit complex is denoted by /3/, for 1 < ? < 3 (see[2l). 

Since each subunit of chaperonin might be empty (denoted by 0) or ATP bound (denoted by 1), 
for the double-ring chaperonin TRiC/CCT with M = 8 + 8 subunits (generally, M is between 14 
and 18), there are altogether 2^ states, denoted hy Sk,k= 1,2, ■ ■ ■2^[1|. The total potential of the 
system with chaperonin in state s^, including potential of the chaperonin subunits and potential of 
the A'^ ATP molecules which are bound to or are around enzyme chaperonin, is given as follows, 

+ 7oln(A^-iVo-M-A^2) (1) 

+ Yi{N-Nq-Ni-N2) 
+ £{M-No-Ni-N2), 

where Ni = N^ +Ni ' for z = 0, 1, 2, A^^ '' (or A^^ is the number of ATP bound subunits which 
have / adjacent ATP bound subunits (in the same ring), and its mirror subunit is empty (or ATP 
bound). 

Based on Boltzmann's law, the probability of chaperonin in state Sk is 



exp{-Ek/kBT) 
YJ^iexp{-Ei/kBT) 



Pk = ^::^ , ^ „ ^. - (2) 



'For any 1 < A; < 2^, si^ denotes the state {/1/2 ■ ■ -'m} of chaperonin which satisfies Ylf=i h^^ ^ —k— I, where 
ij = or 1 for any I < j <M. 



(3) 



Since the potential E^ can be reformulated as follows 

Ek ={ao-£-ri)N^o^ + (ai -e -ri)A^[^^ 
+ (a2 - e - 7i)A^2^^ + (/3i - e - ri)A^o^^ 
+ (j82-e-7i)M'^ + (/33-e-ri)A^f 
+ 7o ln(A^ - A^o - A^i - A^2) + YiN + eM, 

for a given system (i.e. with fixed M and A'^), the expression ?? of probability p]^ can be simplified 
as follows 



exp{-Ek/kBT) 
Lti^M-Ei/kBT) 






where 



and 



+ ^i<)+M^^)+^3A^f (5) 

d,- = a, - £ - 7i, A- = /3, - £ - 7i. (6) 



One can easily show that, the probability Pn that one single chaperonin is bound by n ATP molecules 
can be obtained as follows 

Pn^Y^Pk^ n = 0,l,---M, (7) 

fce/„ 



where 

( M 



£Z, = n, Z, = Oorl \. (8) 



M 

, /=i ('=1 

Obviously, there are ( ) numbers in set /„, so intuitively, the maximum of probability P„ will be 
attained at a number around M/2 [specially, if a, = j8,- = 70 = 0, Pn is binomial distributed, i.e. 
Pn = 2^^( ) . For further details, see the Langmuir adsorption model-^^^]. The average number 



of ATP molecules bound to one Chaperonin is 

M 2^ 

(n)=£nP„=£n(;t);?,, (9) 

where 

M M 

n{k) = Y,h, foryt=£//2'-i + l. (10) 

With the model parameters listed m\T\ theoretical predications of average value (n) and proba- 
bility Pfi are plotted in|3]and|4l which imply that this simple statistical thermodynamics model can 
fit the experimental data of chaperonin TRiC/CCT well.- If chaperonin TRiC/CCT is incubated 
in low [ATP] solutions, the average value (n) increases rapidly with [ATP], but in saturating ATP 
solutions it tends to a constant 8. For convenience, we denote the state of chaperonin TRiC/CCT 
by [«i?2 • ■ ■ hh] [i[i2 ■ ■ ■ i'-ji'sl^ and define Pn* = maxo<n<i6Pn, P,t* = maxi<„<i6/'„. The plots in [5] 
indicate that, for low ATP concentration, the dominant state is [0 • ■ ■ 0] [0 ■ ■ ■ 0] , but for intermediate 
and high ATP concentration ([ATP]> 1 18/iM) the chaperonin TRiC/CCT most likely binds n = 8 
ATP molecules. By further numerical calculations, we found that for high [ATP], the dominant 
state of TRiC/CCT is [11111111] [00000000] or [00000000] [11111111]. One also can see from 
the figures that, although in saturating ATP solutions the average value (n) tends to 8 (see [3]), the 
probability Pn=g tends to a limit which is less than 0.4 (see |4] and [5]). 

Since l&i < 2/3i < 26cq (see[l]), if there is only one subunit of the present TRiC/CCT which is 
ATP bound, the new arrival ATP molecule will most likely attach to one of its adjacent subunits 
(in the same ring), or for simplicity, denote such process as follows [010- • •0][000- ■ -0] + 1 — > 
[01 1 • • ■ 0] [000 ■ ■ ■ 0] or [1 10 ■ ■ ■ 0] [000 ■ ■ • 0] . Similarly, one can easily show that, [01 100000] [000 • ■ ■ 0] + 
1 ^ [01 110000] [000 ■■■0] or [1 1100000] [000 ■ --O] (since 2ai + d2 <j8i+j82 + ai < 2ai + Oq), 
and [01110000][000---0] + l ^ [01 11 1000] [000 ■■•0] or [1 11 10000] [000 ■■ ■0](since 2ai +2d2 < 
«! + 052 + j8i + j32 < 2ai + j8i + J83), etc. One can also verify that for k>4, l&i + {k-2)a2< 
dci + {k-3)a2+$i+$2< 2ai + (A:-4)o;2 + j8i +$3. Therefore, if bound ATP number 2 < n < 8, 
the most likely state of TRiC/CCT is that, all these n ATP molecules are in the same ring and 



adjacent to each other. Moreover, one can easily show that for any 2 < « < 16, all the bound 
ATP molecules will most likely gather together (see|6]). Since 6:2 < ai < 6cq, the ATP binding to 
the same ring is positively cooperated. However, &2 < 0, 7o < but j8/ > implies that the ATP 
binding to different rings is negatively cooperated. 

Conclusions 

In conclusion, from the viewpoint of statistical thermodynamics, a simple method to describing 
the ATP binding process to multisubunit enzymes is present. Each subunit of enzymes might be in 
empty or ATP bound state, so there are altogether 2^ states for one single enzyme with M subunits. 
The potential of each subunit depends not only on ATP binding to itself but also on ATP binding to 
its adjacent (and mirror) subunits. The probability of enzyme in a special state can be obtained by 
the Boltzmann's distribution. Consequently, the probability of bound ATP number in one single 
enzyme and the corresponding average value of ATP number can be obtained. This model can 
fit the recent experimental data for chaperonin TRiC/CCT well.- We find that, in saturating ATP 
solutions, the average number of ATP bound to one single TRiC/CCT tends to 8, and the maximum 
of probability P,, that there are n ATP bound to one single TRiC/CCT tends to a value less than 0.4. 
The most likely state of TRiC/CCT is that, all the bound ATP molecules are adjacent (or mirror) 
to each other, and for ATP number 2 < n < 8, all the bound ATP molecules will most likely be in 
the same ring. The cooperation of ATP binding to the same ring of TRiC/CCT is positive, but it is 
negative for ATP hinging to the different rings. 
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Figure 1: Schematic depiction of the chaperonin TRiC/CCT. TRiC consists of two stacked rings, 
each is made of eight CCT subunits, simply denoted by 1,2, ■■ -8 here. Each CCT subunit might 
be in two states, (empty) or 1 (ATP bound). So there are altogether 2^^ = 65536 states. For 
convenience, the TRiC can also be regarded as two stacked chains with same period 8. The two 
CCT subunits with the same label number are said to be mirror subunits to each other. 
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Figure 2: Potential of one ATP-subunit complex. We assume that, the potential changes of one 
ATP-subunit complex are the same when ATP binds to its left and right adjacent subunits [see (b) 
(b') and (e) (e')]. But they are different from the corresponding change when ATP binds to the 
mirror subunit [see (b) (d) and (c) (e): cui 7^ /3i, a2 7^ J^i]- The a,,/3,- below each figures are only 
potentials of the up-middle ATP-subunit complex, not the total potential of the four subunits. 
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Figure 3: The statistical thermodynamics model prediction [solid curve, see ??], and experimental 
data (solid dots) measured by Jiang et al in,^ of the average ATP number (n) bound to one single 
chaperonin TRiC/CCT. In dilute ATP solutions, (n) increases rapidly with [ATP], and in saturating 
ATP solutions, {n) tends to 8. See[T]for the parameter values used in the theoretical calculations. 
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Figure 4: The model prediction [narrow bar, see ??] and experimental data (wide bar, measured by 
Jiang et al in^) of the probability Pn of bound ATP number in one single chaperonin TRiC/CCT. 
For high [ATP], ?„ gets its maximum at « = 8. See[I]for the parameter values used in the theoretical 
calculations. 
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Figure 5: The maximal values P„* and Pn* of probability Pn of chaperonin TRiC/CCT, and the 
corresponding ATP binding number Hq and/7*. Where P„* =maxo<„<i6^n andP„* =maxi<„<i6/',j. 
The critical [ATP] value in the left figure is about 107 /iM. See [U for the parameter values used in 
the theoretical calculations. 
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Figure 6: Schematic depiction of the most likely configuration of chaperonin TRiC/CCT if it has 
n ATP bound subunits, for n = 2, 8, 10, 15. But one should keep in mind that, if the chaperonin is 
incubated in high [ATP] solutions, it will most likely be bound with n = S ATP molecules [see |4] 
and [51. 
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Table 1: Model parameters used in the theoretical plots in |3l |4] and [51 In the calculations, we use 
N = No = 25 for the chaperonin TRiC/CCT incubated in 5/iM ATP solution, and A^ = A^o[ATP]/5 
for other ATP solutions. All the following parameter values and A'o = 25 are obtained by fitting ?? 
?? to the experimental data obtained by Jiang et al.'^ The unit of parameters d,-, j8/ is hgT. 
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